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We consider a section of a half-filled chain of free electrons and its entanglement 
with the rest of the system in the presence of one or two interface defects. We find a 
logarithmic behaviour of the entanglement entropy with constants depending continuously 
on the defect strength. 



The entanglement of different parts of a quantum system in the total wave function can be viewed as 
the result of a coupling across their common interface. Correspondingly, the entanglement entropy turns 
out to be proportional to the interface area in simple models [1-3]. Modifications of the interface will 
therefore modify the entanglement and in one-dimensional systems with short-range interactions, where 
the interface reduces to a point, a single defect is expected to have a marked influence. This was recently 
pointed out by Levine [4] who applied this idea to the case of a Luttinger liquid with one impurity. 
He obtained within perturbation theory a reduction of the entanglement entropy S which is strong for 
repulsive and weak for attractive interactions, if one considers a large subsystem. Roughly speaking, this 
corresponds to the known influence of interactions on the impurity strength in this case [5]. Levine's 
approach using bosonization is interesting, but it also suggests to look directly at the simplest case and 
to investigate the problem for non-interacting electrons with specific defects. 

In the present paper we therefore study free electrons hopping on a chain for the case of half filling. In 
spin language, this corresponds to an XX model. This is a critical system and the entanglement entropy 
between a subsystem of length L and the rest is given by the conformal result [6-10] 

S= C -\nL + k (1) 

with the value c = 1 for the central charge. How does this change in the presence of a defect at the chosen 
interface ? For one case, the answer is known : if the defect is such that it cuts the chain, the subsystem 
has a free end and thus only one connection to the rest remains. This changes c — > c e // = c/2 and 
k — > k' , but the logarithm is unaffected [7]. In the following we show that the same holds for arbitrary 
defects, either at one or at both interfaces. These defects can be changed bonds or changed site energies. 
The calculations are numerical and based on the determination of the reduced density matrix p from the 
one-particle correlation function of the total chain [11,8,12]. The necessary formulae are given in Section 
I. Eigenvalue spectra of p for various defects, which differ in a characteristic way from the spectrum of 
the homogeneous system, are presented in Section II. From the eigenvalues, the entanglement entropy 
S = —tr(plnp) is obtained and discussed in Section III. We determine the effective central charge c e // 
for bond and site defects and discuss its behaviour in terms of a simple model. In contrast to the constant 
k, it is always reduced by a defect and therefore the entanglement always becomes smaller if the size L 
of the subsystem is large enough. A brief summary is given in Section IV. 

I. BASIC FORMULAE 

We consider a system of free fermions hopping between neighbouring sites of an infinite linear chain. 
The corresponding Hamiltonian reads 
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H = -^i„(4c«+l +4 +1 C„) +J2 A ncU n (2) 

n n 

where t n is the hopping matrix element and the 'hat' denotes quantities of the total system. In the 
following, we set t n = 1 and A„ = except at the boundaries of the subsystem which consists of the 
sites i = 1, 2, ...L. We will mainly consider the case of one bond defect at the left boundary, t — t, but 
we will also give results for two equal bond defects at the two boundaries, to = th = t, and one or two 
site defects next to the boundary, Ai = A or Ai = Ai = A. 

The total system is assumed to be half filled and in its ground state |0 >. The reduced density matrix 
then has the form [11,14] 

p = K exp (- H H c \ c i) ( 3 ) 

where K, is a normalization constant and the matrix Hij follows from the one-particle correlation function 
of the total system 

Cmn =< 0| <4c„ |0 > (4) 

via the relation 

H = ln[(l -C)/C] (5) 

Here C denotes the L x L submatrix of C with the sites restricted to the subsystem. For a homogeneous 
infinite system, C is given by 

With defects, the translational invariance is lost and for single defects C has the general form 

C l3 =C°(i- J )-C 1 (i+ J ) (7) 

Thus Cij is the difference of a Toeplitz matrix depending on (i — j) and a Hankel matrix depending on 
(i + j). Physically, the term C 1 leads to oscillatory behaviour of the correlations as one moves along the 
chain. For site defects, these are the Ruderman- Kittel oscillations if one looks at the density Cu. For 
bond defects, the density is unaffected but the oscillations appear e.g. in the nearest-neighbour (bond) 
correlation function Cis + \. The quantity C 1 can be obtained by using the scattering phase shifts and 
including the contribution of possible localized states caused by the defect. Thus a single weak bond 
t = e v < 1 in an infinite total system leads to 

C\l) = -isinh^ (e-"J, - e v h_ 2 ) (8) 



where 



Jo 71 smh i/ + sm q 

For t = 0, this simplifies to C 1 (/) = C°(l) and C assumes the known form for a system with an open end. 
For a strong bond, t > 1, one has to add the contribution 

C 1 ^) = -sinhzy exp(-v(l - 1)) (10) 
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from the occupied bound state at energy —2 cosh v below the band. Similarly, a site defect A = 2 sinh v > 
leads to 



C l (l) = (/,_! - V 3 + 2smh^/ i _ 2 ) 



(11) 



By diagonalizing C one obtains its eigenvalues Cfe (0 < Cfc < 1)> from which the eigenvalues s k = 
ln((l — CO/CO 01 ^ (— oo < e k < oo ) follow. The reduced density matrix then takes the diagonal form 



p = /Ccxp (- ^2e k c\c k ) 



(12) 



The entanglement entropy is defined by S = —tr{p\n p) (we are using the natural logarithm, not the one 
to basis 2) and reads in terms of the ( k 



S = - £ [Cfc InCfc + (1 - Cfc) Mi - CO] 



(13) 



and in terms of the e k 



In [1 + exp(-Sh)] 



exp(e fe ) + 1 



(14) 



Only the ( k which are not too close to or 1 rcsp. the e k which are not much larger than 1 in magnitude 
give a sizeable contribution to the entropy. 



II. DENSITY-MATRIX SPECTRA 



We have calculated the correlation matrix C for the case of single defects by evaluating the integrals 
for C 1 numerically. From this the low-lying single-particle eigenvalues e k with \e k \ < 25 were obtained. 
Higher values cannot be reached with standard double precision routines because the corresponding Q k 
lie too close to or 1. However, as noted above, they are unimportant for S. In the following we always 
consider even L. 

Spectra for one bond defect are shown in Fig. 1 for a subsystem of L — 50 sites. The eigenvalues come 
in pairs (e, — e) due to the particle-hole symmetry of the problem. For a homogeneous system, t = 1, one 
finds the slightly bent curve known from previous investigations [15,13]. 

For a weak defect, one can see a shift of the dispersion curve which is upward for positive e k and 
downward for negative ones. This shift becomes stronger as the bond becomes weaker. In addition, 
oscillations appear which also increase for weaker bonds. Most importantly, however, more and more of 
the low £fc follow a steeper dispersion curve as t approaches zero. This curve has a slope approximately 
twice as large as for the homogeneous system and no wiggles and represents the result one obtains directly 
for the system with a free end. Thus the spectrum with the defect is a mixture of those for the two cases 
t = 1 and t = 0. 

For a strong defect, the situation is basically the same. However, here the two lowest eigenvalues e 
and E\ play a special role because they approach zero as t — > oo . This moves the left and right part of 
the spectra two units apart. This feature can be understood from the correlations in this limit. The two 
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FIG. 1. Low- lying single-particle eigenvalues for a subsystem of L=50 sites with one boundary defect. 
Left : Weak defects, t = 1; 1CP 1 ; 1CP 2 ; 1CP 3 ; 1CP 4 (from bottom to top in the right part of the figure). 
Right : Strong defects, t = 10; 10 2 ; 10 3 ; 10 4 (also from bottom to top). The lines are guides for the eye. 
The numbering is such that positive e k have positive values of k. 

localized states below and above the band then exhaust the local Hilbert space at the two bond sites 
and effectively cut the system. Thus C\j = for j > 1 and the diagonal term C\\ = 1/2 gives one 
eigenvalue zero. The other one appears because the remaining part of the subsystem has an odd number 
(L — 1) of sites. It is absent, if the full subsystem has L + 1 sites. Then the spectrum for i ^ 1 is exactly 
the same as for L sites and 1/t <C 1 up to the remaining zero eigenvalue. 




FIG. 2. Low-lying single-particle eigenvalues for a subsystem of L — 50 sites with two boundary defects. 
Left : weak bond defects. Right : strong bond defects. The parameters are the same as in Fig. 1. 
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Spectra for two bond defects are shown in Fig. 2. These were calculated by finding the eigenfunctions 
of H for a total system of size N = 1000 numerically. This introduces small finite-size effects but these 
are not visible on the scale of the figure. There is a shift of the curves as for a single defect, but now it 
affects all eigenvalues if the defects are weak. Thus a gap is opened in the spectrum. For t — > all Sk 
diverge. This is reflects the fact that the subsystem is, in this limit, completely decoupled from the rest 
and the ground state |0 > becomes a product state. Then p has one eigenvalue 1 if all single-particle 
levels arc unoccupied and all other eigenvalues vanish. If the defects are strong, two Sk go to zero while 
the remaining eigenvalues again diverge. The zero eigenvalues are related to the localized states at the 
two bonds which belong to the subsystem as well as to the environment. The ground state |0 > then has 
four components which differ by the location of the electrons on the two bonds and the four resulting 
eigenvalues 1/4 of p are produced by the two vanishing 

For comparison, the spectra for two site defects are shown in Fig. 2. These defects destroy the particle- 
hole symmetry of the problem and therefore the spectrum of the Sk is in general no longer symmetric 
with respect to zero. Apart from that, it shows the same features as for two bond defects. In particular, 
also the site defects cut off the subsystem if the site energy A becomes large and thus all except two 
diverge for A — > oo. In this limit, the reflection symmetry of the spectrum is also restored. 
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FIG. 3. Low-lying single-particle eigenvalues for a subsystem of L = 50 sites with two site defects at the 
boundaries. The site energies are A = 0; 2; 5; 10; 100 from bottom to top in the right part of the figure. 

III. ENTANGLEMENT ENTROPY 

Using the values of the or the £fc, one can obtain the the entanglement entropy from (13), (14). 
Results for single bond defects are shown in Fig. 4 for system sizes between L = 20 and L = 100. 

The plot on the right hand side shows that S varies logarithmically with L in all cases. There is no 
indication of a In 2 L-tcrm as found in [4] . There is only a small variation of the slope with L and its 
asymptotic value can be determined very well by an extrapolation in 1/L. It is remarkable that one can 
see this logarithmic behaviour so well, since for all sizes considered here only about 20 single-particle 
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eigenvalues contribute if the accuracy is set to 1CP 10 . Thus one is far away from a limit in which the 
eigenvalue spectrum becomes dense and where the logarithmic dependence could result from the density 
of states. 
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FIG. 4. Entanglement entropy as a function of L for one bond defect. From bottom to top: t=0.1; 0.5; 1.0; 2.0 

Thus we find that the entanglement entropy has the same form (1) as for the homogeneous system, 
but with an effective value c e // which depends on the strength of the defect. The same holds for two 
bond defects and also for site defects (see below). We first discuss the case of one bond defect further. 

The asymptotic value of c e ff, determined by extrapolating the data from L = 20 — 100, is shown in 
Fig. 5 as a function of the defect strength. 
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FIG. 5. Effective central charge as a function of the bond strength for one bond defect. Full line : numerical 
result, dashed line : approximation Eqn. (16). 
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One sees that it varies smoothly between the limits c e // = 1/2 for a subsystem with a free end and 
c e ff = 1 for a homogeneous system. The figure shows c e // only for t < 1 because one finds that c e // is 
symmetric under t — > 1/t. 

One can try to find a simple analytic function which fits the numerical curve in Fig. 5. Since the defect 
is characterized by the scattering phases, one could expect that c e // depends on these. Because one is 
dealing with an asymptotic property the phase shifts at the Fermi surface should enter. These are given 

by 

<5± = ±[2arctg(t)-|] (15) 

for functions symmetric and antisymmetric about the defect, respectively. One is then lead to the 
expression 

c e// = \ + \[~ arctg(i) - arctg(l A) ] 2 (16) 

The expression in the brackets also appears (as an exponent) if one calculates the time autocorrelation 
function of a spin in an XX spin chain with a bond defect [16]. The function (16) has the correct limits 
for t = and t — 1 and is symmetric under t — > 1/t. As can be seen from the figure, it approximates the 
numerical data quite well. The deviations are less than 10~ 2 , but the agreement is not perfect. Especially 
for t — > 0, the data do not seem to vary quadratically in t but with a power 1.8. Replacing the power 2 
in (16) by this value improves the fit over the whole intervall substantially. The deviations become less 
than 10~ 3 and are almost invisible on the scale of the figure. 



The two terms in (16) can be viewed as the contributions of the unperturbed and the perturbed 
interface, respectively. One can understand this result from the form of the single-particle spectrum. In 
the previous section it was found that it consists of two pieces. If one assumes that in the limit L — > oo 
the levels in these two branches become equidistant and dense such that (for positive Sk) 



2k a/ InL, e k < a 
k a/ InL, Ek > a 



(17) 



where a = a(t) depends on t such that a = for t — 1 and a = oo for t = and a is a constant, the 
entanglement entropy can be written as 



S = 2\nL 



1 

2a 



de s(e) + 



de s(e) 



(18) 



Here s(e) is given by the bracket in (14) and the factor of 2 in front takes care of the negative levels. 
Denoting the second integral by G(a), one obtains 



£ = 2 InL [I(G(0)-G(a))+G(a)] =InL [G(0) + G(a)] 



(19) 



Thus each interface brings its own amplitude for the logarithm : the unperturbed one G(0) and the 
perturbed one G(a). For the homogeneous system the bracket must have the value 1/3, from which 
G(0) = 1/6 follows. Therefore (19) gives for c e // 



c e ff = i + 3 G{a) 



(20) 



If the exact dependence of a on t were known, one could evaluate the integral G(a) and determine c e ff 
in this way. The oscillations of the Sk which are neglected in (17) could also be included. 
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The discussion so far has dealt with one defect. However, one can apply it also to the case of two defects. 
The results of Section III show that then the lowest eigenvalues are just missing and one also finds that 
the gap in the spectrum for two defects corresponds roughly to the bending point in the spectrum for 
one defect. Therefore one has only the second region in (17) and finds instead of (19) 



S = lnL[2 G(a)} 

from which the value of c e / / can be read off. This also gives the relation 



„( 2 ) 

-eff 



2c (1) -1 



(21) 



(22) 



between the effective c- values for two defects and for one defect. Due to the lower accuracy of the calcu- 
lations for two defects, this relation could not be checked precisely, but the deviations are less than some 
percent and vanish near t = and t = 1. 

Results for the case of one site defect are shown in Fig. 6. Here the entanglement entropy is symmetric 
under A — > — A and only positive values of A have to be considered. One sees that c e // decreases again 
from the value 1 to 1/2 as the perturbation becomes stronger, because the defect cuts the chain in the 
limit A — > oo and thus acts like a bond with t = or t = oo. 
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FIG. 6. Effective central charge as a function of the site energy A for one site defect. The values were obtained 
by extrapolation in 1/L 



Finally, we present results for the constant k in the entropy, see (1). They were obtained by extrap- 
olating S(L) — (c e ///3) InL from the values L = 20 — 100. As Fig. 7 shows, this quantity k is not 
symmetric under t — > 1/t for one bond defect. Rather it increases from about 0.5 for t = to about 
1.2 for large values of t, with a small initial dip (which depends sensitively on c e //). At t = 1 we find a 
value k = 0.726 which is in agreement with the findings in [9] but differs slightly from the value obtained 
from the asymptotic analysis in [10]. The limiting value for t — > oo is found to differ from the value at 
t = by In 2. If one considers systems with an odd number of sites L + 1 for t 3> 1, this also holds 
and follows in this case directly from the identity of the single-particle spectra up to one zero eigenvalue 



8 



(which contributes In 2) noted in Section III. The relatively large value of k can make the entropy with a 
defect larger than without one if the subsystem is small enough (compare Fig. 4). However, in the limit 
L — > oo the effect of the reduced c e // always dominates and leads to a reduction of the entanglement. 
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FIG. 7. Constant k in the entanglement entropy as a function of the bond strength t for one bond defect. The 
values were obtained by extrapolation in 1/L 



We have investigated a chain of free electrons in its ground state at half filling. The entanglement 
properties were calculated for the case of one or two defects located at the contact points between the 
subsystem and the rest of the chain. The density-matrix spectra were seen to differ in a characteristic 
way from those for a homogeneous system. Nevertheless, the logarithmic dependence of the entanglement 
entropy on the size persists, only the amplitude changes and becomes dependent on the defect strength. 
The mechanism could be described by a simple model based on the features of the single-particle spectra. 
A simple approximate formula for the effective central charge was also given. Although we focussed on 
bond defects, the same features are found for site defects. The logarithmic law means that the system 
remains critical, as expected for a localized perturbation. In the two-dimensional system associated with 
the quantum chain, the defect becomes a line as noted and used in [4]. Such a two-dimensional system 
with a cut (because one considers the reduced density matrix) and one or two straight defect lines per- 
pendicular to this cut can be mapped conformally to a strip [14] but the defect lines are bent by the 
mapping. For a single defect, the line runs only in the left or right half of the strip, coming from one point 
on the boundary and going back to another one. For two defects, the two lines are attached to different 
edges and touch in the center of the strip. The more complicated density-matrix spectra found here must 
be related to this feature. It would be interesting to derive them analytically. Also one wonders if the 
asymptotic formulae used in [10,17] to calculate S can be generalized to the present case. 
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IV. CONCLUSION 
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